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Abstract 



Cusps in superconducting cosmic strings produce strongly beamed electro- 
magnetic radiation. To calculate the energy emitted requires taking into 
account the effect of the charge carriers on the string motion, which has pre- 
viously been done only heuristically. Here, we use the known exact solution 
to the equations of motion for the case where the current is chiral to update 
previous calculations for the total energy, spectrum and angular distribution 
in that case. We analyze the dependence of the radiated energy on the cusp 
parameters, and discuss which types of cusp dominate the total radiation 
emitted from an ensemble. 
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I. INTRODUCTION 



Cosmic strings are one dimensional topological defects which may have been left behind 
by phase transitions in the early universe |I]J2]. in the zero width approximation we can 
describe the dynamics of the strings by the Nambu-Goto (NG) action H|J, namely 



/ d 2 (V^ (1) 



where /i is the energy per unit length of the string and 7 the determinant of the worldsheet 
metric parameterized by the coordinates (° and ( 1 . Due to the reparameterization invariance 
of the Nambu-Goto action, we have the freedom to impose several conditions on the functions 
that denote the string position, in other words to fix the gauge. The most convenient gauge 
for ordinary strings in flat spacetime is the conformal gauge, in which we impose 

|xV,r)| 2 + |x(a,r)| 2 = l (2a) 
x'(cr, r) • x((7, r) = (2b) 

where we have fixed (° = t = x° and C 1 = cr, and primes and dots denote differentiation with 
respect to a and r respectively. . With these conditions the equations of motion become 

x>,r) =x(a,r). (3) 

The general solution is 

x(<7,T) = -(a(<7-T)+b(<7 + r)) (4) 

where a and b are arbitrary functions constrained by the gauge conditions to satisfy |a'| = 
|b'| = 1. From Eq. (|j) we see that 

x'=i[aV-r) + bV + r)] (5) 

and 

1 

x = - 
2 



a'( 



a - t) + b'(<7 + t)) . (6) 



If there is a point along the string at which a' = — b', then x' = and even more importantly 
|x| = 1. This means that there are points on the string that move at the speed of light and 
where the string doubles back on itself. (See Fig. |l|.) These events are called cusps 0. 

It is clear that this kind of dynamics falls outside of the domain of validity of the as- 
sumptions involved to obtain the Nambu-Goto action. Nevertheless it has been shown by 
field theory simulations that the motion of the string is quite accurately described by the 
NG solution up to the point of interaction between the two branches of the string ||. 

The violent motion of the string at the cusps suggests the possibility of a range of effects 
which could have prospects of detection ||7|-|T5|j. 

In 1985 Witten |16j showed that certain particle physics models could lead to the forma- 



tion of strings with superconducting behavior. It was soon realized that these models could 
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FIG. 1. Form of the cusp region for ordinary and superconducting cosmic strings: pure 
Nambu-Goto cusp (solid); including charge carrier back-reaction (dashed). 



in principle have a much richer phenomenology due to the inertia of the charge carriers and 
their coupling to the long-range electromagnetic field. 

In particular we can think that if the current on the string is not very large, then the 
string motion would not be affected, and calculate the electromagnetic radiation from a 
Nambu-Goto prescribed motion for the string. However, the power emitted by radiation 
near the cusp in this approximation is divergent due to the infinite Lorentz factor of the 
NG cusp P,[T7| . This infinite result suggests that electromagnetic back-reaction is crucial in 
order to understand the superconducting string cusp. 

The short-range effects of the electromagnetic field can be taken into account by renor- 
malization of the charge carrier inertia and charge |!6| . |l8| f| . The charge carrier inertia term 
prevents the appearance of an infinite Lorentz factor, and thus the calculation yields a finite 
result. In principle, this result could be further modified by back-reaction due to long-range 
fields (i.e. radiation), but for generic values of the parameters we expect this correction to 
be small |19|]. Further approximating that the current is chiral, at least in the vicinity of 
the cusp, we are able to do a proper calculation of the total power radiated. 



II. SUPERCONDUCTING STRINGS 



The equations of motion for a superconducting string can be written |16 

d a (V=7 (^ ab + ab ) < 6 ) = - V=7 

d a (V=TT a V,6) = -\qe ab F^ a xl 
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(7a) 

(7b) 
(7c) 
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where /i is the energy per unit length of the string, q the electromagnetic coupling of the 
fields on the string, 6 ab the energy-momentum tensor of the charge carriers, 

e ab = 7 ac 7 M 0,c <f>,d - \ l a V d <t>,e (t>,d , (8) 
the auxiliary scalar field in terms of which we can express the conserved worldsheet current 

J a = g-^=e a6 0, 6 , (9) 

and j M the 4-dimensional current, 

f( x ) = [ d 2 c^ r x^ 4 \x - x(0) ■ (io) 



These equations of motion are obtained from a two-dimensional effective action for the 
superconducting string and they are valid as long the current, measured in an invariant way, 
is not too large, 

JaJ a ^\Jcrlt\ 2 , (11) 

where J cr i t is the limiting current at which charge carriers can be ejected. 

The leading effect of electromagnetic self-interaction of the current-carrying string can 
be taken into account by a redefinition of the auxiliary field (f> and its charge ]TB| , P^| ,|2| . This 



enables us to include the near-field effect on the left-hand side of the equations of motion. 
If we take this as our only back- reaction effect, then we can rewrite the equations above in 
a much simpler way by taking q and <fi to be the renormalized quantities, to get 

d a (V^l (m ab + 9 ab ) aft) = (12a) 

d a (V=T7°V, 6 ) = (12b) 
d 1/ F tiV = -Anf (12c) 

This procedure decouples the first two equations from the third one, and they now give 
the motion of a 'neutral' superconducting string. In the last equation, F^ u is the external 
field emitted by the current moving according to the first two equations. 

Numerical simulations have suggested that as a superconducting cosmic string 

loop contracts under its own tension, it would approach the 'chiral' limit, in which the 
current on the string is null, in other words, J a J a = 0. This, in turn, means that the 
auxiliary scalar field fulfills the condition, 

7 a V,60,a = O. (13) 



Furthermore, it can be shown |18[ that if the current is not chiral it would become 
supercritical around the cusp and charge carriers would be ejected. We therefore expect 
that this process would drive this region of string around the cusp to the chiral limit. These 
arguments suggest that most of the electromagnetic radiation would come from chiral current 
carrying string cusps. 
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If the current is chiral, the set of equations describing the string and the 'neutral' current 



can be solved completely [P2,1S, P31 . In this case, the general solution for the scalar field is 

(f>(a,T)=F(a + T), (14) 
where F is an arbitrary function, and the string position is given by 

x° = r (15a) 
x=i[a(<7-r) + b(<7 + T)], (15b) 

with constraints for the otherwise arbitrary functions a' and b', 

|a'| 2 = 1 (16a) 
4 F' 2 

|b'| 2 = l-^— . (16b) 

This description of the chiral strings makes evident that they behave very similarly to 
the Nambu-Goto case. As before, the functions a' and b' live on the surface of spheres (if the 
current is constant), but now their radii are different for the two cases. This change, however, 
brings important consequences for the string motion. As the magnitude of the current 
increases, the average velocity of the loop decreases, arriving at the limiting stationary case 
when 4F' 2 //i = 1, the chiral vorton |20|j . 

This also means that the string can only reach a finite Lorentz factor at the cusp-like 
regions, in contrast with the Nambu-Goto case. This effect can be easily understood in terms 
of the resistance of the charge carriers to being accelerated by the string tension. Equations 
(|T5j) also imply that |x'| > at the cusp so the actual shape of the string around the cusp 
is rounded off by the effect of the charge carriers, as shown in Fig. |I[ (Another way for the 
string to be rounded off, producing a self intersection, is discussed in |^4| .) 



III. ELECTROMAGNETIC RADIATION FROM THE CUSP 

Using the solution for the string motion from the previous section, which as discussed 
is free from divergences in its Lorentz factor, we can now compute, using Eq. ( |12c| ), the 
electromagnetic radiation from the cusp-like region. Throughout this section we follow 



closely the calculation of the radiation by Spergel et al. iT7| , and introduce the necessary 
modifications for the chiral string motion. 

Since the radiation is mainly produced by the region with a large Lorentz factor, we can 
estimate the energy output from the position of the string expanded in a Taylor series around 
the cusp, which we take as a = and r = 0. We also assume throughout the calculation 
that 2 \F'\/ ^JJi <C 1, since otherwise the string would not reach high Lorentz factors at all. 

Using the light-like variables cr+ = a + r and o~_ = a — r, we can expand the functions 
a and b near a± = as 

a(<7_) = a[,o-_ + l&'y_ + -4'a 3 _ + ■■■ (17a) 
2 6 

b(<r + ) = b' a + + \h'y + + ibgV» + ■ • • , (17b) 
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where we center our coordinate system at the position of the cusp, so that x = 0. In 
analogy with the Nambu-Goto case, we consider the cusp to be the place at which a' and b' 
are anti-parallel, 



b' 



(18) 



which corresponds to the point of maximum concentration of energy and the minimum value 
of |x'| ||19|| . Assuming that the current remains constant over the region of interest, we can 
see that the constraint equations fll6|) impose a constant magnitude for the vectors a' and 
b'. We can choose for simplicity 



b' 



—x 
Ax 



(19a) 
(19b) 



where A = 1 — e, with e < 1 and constant. Due to the fact that |a'| and |b'| are constants 
we have the following relations, 



a • s!q 

/ w 
a ■ a 

K ■ K 



b'o-b^O 



\K\ 



(20a) 
(20b) 
(20c) 



Using Eqs. (plj, |2(|) we can express the most general expansion of the string position near 
the cusp up to third order in a± in the form 



a <7_ 



bfo-. 



a 2 a 3 



6 

A(T 4 



— (7_J X 

(l 2 + 7 2 ) 
6A 



1 2 1 3 

- a a H — no 

2 ~ 6 ^ ~ 



5 a 3 



(21a) 



x 



1 



/3ai + -wai) y+ (- 7 < + ^<) z. (21b) 
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The parameters a, j3, 7 are of order L _1 , with L denoting the typical length scale of the loop, 
and 6, zu, g, r\ are of order L~ 2 . 

Using standard results in electromagnetic theory [^5], we can calculate the energy radi- 
ated per unit frequency and per solid angle for a source given by Eq. ( |TUD using the expression 



17 



d 2 E r , 



rad 



U) 



J da J dr{n x c) exp[— iu;(ri'*x /4 )] 



du dVt 4tt 2 

where is the unit null vector from the source point to the observation point, 

e 2 



1, cos 9, sin 6 cos ip, sin 6 sin (fj ~ ^1, 1 — — , 9 cos </?, 6 sin (f^j . 



(22) 



(23) 



This last approximation is justified since the radiation is highly beamed in the direction of 
movement of the cusp, in this case, x. The four- vector &■ is the source current, which in our 
case takes the form 
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d i = q f ? b d a <f>d b x ft (24) 
where a, b = 1, 2 stand for a and r. In the chiral case, the four-vector current becomes 



c" = q 



' <9x M dx^ 



2qF 



, dx^ 
daZ 



da \ da dr 
Keeping only the lowest terms in the angle 9, we get 

" J da_ C(«r_, 9, if, F') exp[iuj(Ja_ + Ka 2 _ + Lai )} 



d 2 E rad 



did dQ 167T 2 



x 



J da + exp[— %lo (Sa + + Ta 2 ^ + Ucr+)] 



with 



and 



C(er_, 9, ip, F') = qF'i — aO sin ipa^,—9 sin ip, aa_ + 9 cos ip 



K = -a cos 09 
4 ^ 

L - — 
12 



and 



5= o 1 - A ( 1 " V) 



T = — - (/5 cosy? + 7 sin 93) 



U 



12A 



(25) 



(26) 
(27) 

(28a) 
(28b) 
(28c) 

(29a) 

(29b) 
(29c) 



The two integrals in Eq. (^6|) can be evaluated, first changing variables to get rid of 
the term quadratic in a± in the exponent, and then using the following forms of the Airy 
integrals: 



and 



dx cos 



dxx sin 



~£ I x H — x c 



— M x H — of 



(30) 



(31) 



where are the modified Bessel functions of order v. With a little bit of algebra we get 
to the result, 
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d 2 E. 



rad 



UJ- 



doo dVt 16n 2 



(qF'f (~#i(£ + )x + )' 



x 



9 2 sm^i-K^X-) +a 



\3\/3 



(32) 



where 



I3SU-T 2 _ /3[2e + # 2 F]\ 



1/2 



V 3C/ 2 \ 
2 / T 2 N 



1 /[2e + # 2 Y] 
X+ " 6* { /? 2 + 7 2 



3 X 1/2 



and where we have defined the dimensionless quantity, 



Y 



((3 simp — jcosipY 



sin 2 ((y9 - y? ) 



(33a) 
(33b) 

(34) 



where (po is the angle between a^' and b' ', so that tany?o = l//3- On the other hand the 
expression for the same type of parameters for the integral in o~_ are much simpler, 



'3JL - K 2 



X- 



3^ 



3L 2 



oo J 



= Vse 

K 2 \ 



sin ip 



a 



3L J 



X- 



6 



uo6 3 



sin 3 if 



(35a) 
(35b) 



It is important to note that the different behavior of the superconducting chiral string cusp 
is encoded in the presence of the parameter e. If we take the limit e = 0, we would recover 
the pure Nambu-Goto cusp. 

Using the above definitions, we can simplify Eq. (R^) to 



d 2 E rad fqF'sm 2 (p\ 2 



duo dQ 



37Tft 



K) 



2 ( 2e + 6 2 Y ^ 

V P 2 + i 2 , 



+ 



We can now integrate over frequencies with the change of variables z = to get, 



dR 



rad 



dQ 
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7T 



9 3 | sin A (p\ 



ft 



1/3 



((3 2 + ~f 2 f 3 



where we have defined the function T as 



T{a) 



dza^z 2 



Ki (az) Ki(z) 

3 3 



Ki(az) K2(z) 



(36) 



(37) 



(38) 



The integral can be done exactly in hypergeometric functions. The function T has asymp- 
totic behavior 



const a — ► 
const /a a — ► oo 



(39) 



(with different constants). 
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10 



A/h 6 \ J 

2 

tt/4 " tt/2 3tt/4 tt 

^0 

FIG. 2. Numerical coefficient .4 for an m = 1, n = 2 Burden loop. 

In order to integrate Eq. (|37l) with respect to the angle 9, we first note that we can extend 
the range of integration to infinity since only the region around the direction of the cusp 
contributes. In order to extract the dependence on the parameter e, we perform another 
change of variables to 9' = so the final expression for the energy emitted at the cusp 

becomes 

E rad = A( a ,(3, 1 )(qF') 2 e- 1 / 2 , (40) 

where A(a, j3, 7) is a constant which can be written in terms of the magnitudes of and 
b>Q and the angle between them, ip , as 

94 la"! 1 / 3 /-27T /-co 

A =^l d *L "r-'iin^ 

(41) 

Figure [2] plots the value of A/L versus ipo for an m = 1, n — 2 Burden |^B[ loop, which has 
I a.Q I = 2n/L and |bg| = 47r/L at the cusp. 

We can understand the form of A as follows. When 6' <^ 1 (i.e., 9 <^ y/e), the argument 
to T goes as 9'~ 3 sin -3 (p ^> 1, and the integrand becomes proportional to 9', which is just the 
geometrical factor from the integration. Thus the radiated power per solid angle is uniform 
for small angles. Now we consider the case with 9' large. When sin(<£> — ip ) is far from zero, 
then 9' ^> 1 makes the argument of T constant. In such directions, the radiated power is 
suppressed by 9'~ 2 , so their contribution is small. If sin(y9 — (p Q ) is small, then the argument 
of T goes to zero, and thus the value to a constant, unless sin<^ is also small. Thus, as long 
as ipo is far from and 7r, the radiation is confined to a beam with radius 9 ~ ^fe. 

However, if ip is near (or, equivalently, near 7r), then it is possible to have sin((p — ip ) 
and sin ip small simultaneously. In that case, there is a range of 9' up to about 1 , in which 



2 + 9' 2 sin% - <p ) 



3/2 



Ibffl 



0/2 



sin Lf 



O(L) 
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large contributions to the integral are still possible. This is the cause of the divergence seen 
in Fig. ||: A is proportional to y^ 1 f° r 1. In this case, the beam can be greatly 

elongated in the plane of the cusp. In the limiting case, = or = 7r, the Burden loop 
has a cusp at all times, which rotates around the string. In cases close to this, we see that 
the cusp is longer lived than usual and thus has the opportunity to beam radiation through 
a wider angle. 

In this regime, Eq. ( f20"| ) is not correct, because we cannot extend the 9 integration to 
infinity. Instead, we should have a cutoff at 9' ~ e -1 / 2 . If we include this effect, then the 
growth of A with ip^ 1 is cut off at ip Q ~ ^fe. 

Now suppose that we have a set of loops with ip distributed evenly between and it. 
To find the contribution due to the loops with ipo near 0, we can integrate 



dtpo E(ip ) ~ / d(f) (j) ~-lne. (42) 

Thus if e is small, loops with and bg nearly parallel or antiparallel dominate the total 
electromagnetic radiation, but the effect is only logarithmic. 

Returning to the generic case, let's now remember that we have defined e to be the 
deviation of the magnitude b' from unity, 

|b'| = A = 1 - e (43) 

so in the chiral string model, for a current much smaller than the energy per unit length of 
the string, 



/ A\F'\ 2 A\F'\ 2 
|tf| = ,/l-U-*l-±-L (44) 



so 

4IF 



/|2 

(45) 



2/i 

It can be shown |L9| that the parameter e determines the maximum Lorentz factor 
reached near the cusp, namely, 

T max = ^e)- 1 ' 2 (46) 

so we can rewrite Eq fl4"0| ) as 

E rad ~q 2 ^Jl \F'\ L~L {qF'f T max (47) 



in agreement with earlier estimates |p|,|17|. We can also write the total energy output in 
terms of the physical current as 

Erad ~ q J v 7 ^ L ■ (48) 
Equations ([|7],[IS|) are valid when <j) is far from and ir. 
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IV. DISCUSSION 



Electromagnetic radiation from superconducting cosmic strings cusps has traditionally 
been considered a distinctive signature of superconducting string models. Recently this 
interest has been revived by as possible connection with gamma ray bursts [|l4t|l5| and 
ultra-high energy cosmic rays. The total energy output from a pure Nambu-Goto string 
cusp (i.e. no back-reaction included) yields an infinite result due to the infinite Lorentz 
factor of the tip of the string, where the charge carriers get concentrated. It is therefore 



necessary to impose some sort of cutoff for this process. Early calculations [p%17| assumed 
the cutoff to be where the energy in the charge carriers was comparable with the energy of 
the string. This procedure gives a linear dependence of the total energy with the maximum 
Lorentz factor. We redo this calculation in the context of a chiral current on the string. 
This is the most interesting case since loops are driven to chirality by charge carrier ejection 
especially near the cusps ITS . 



The full system of equations for a superconducting string is very complicated. Neverthe- 
less in the case of a chiral current they are much simpler and can be solved exactly |22 , |19 ,23 



This allows us to compute the electromagnetic radiation from a chiral-current string tra- 
jectory taking into account in an exact way the back-reaction of the charge carriers on the 
string motion. In order to do this we closely follow the procedure described in Spergel et 



al. [17] and we extract the dependence of the total energy on the maximum Lorentz factor 



that the string has in the pseudo-cusp region. We found the same linear dependence in 
the maximum Lorentz factor as earlier estimates. This suggests that the intuitive picture 
of thinking about the back-reaction of the charge carriers as a process that only becomes 
important when their energy density is comparable to the energy of the string is correct, 
making it possible to extend this result to all types of currents. 

In the case of cusps where the parameters a' ' and b^' are nearly parallel or antiparallel, 
the radiation can be much larger than in the generic case. For small currents, the radiation 
from such cusps dominates in the total flux, but only by a logarithmic factor. 
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